SPECTRAL ANALYSIS OF ONE-DIMENSIONAL HIGH-CONTRAST 
ELLIPTIC PROBLEMS WITH PERIODIC COEFFICIENTS 

K. D. CHEREDNICHENKO* , S. COOPERt, AND S. GUENNEAU* 

Abstract. We study the behaviour of the spectrum of a family of one-dimensional operators with 
periodic high-contrast coefficients as the period goes to zero, which may represent e.g. the elastic 
or electromagnetic response of a two-component composite medium. Compared to the standard 
operators with moderate contrast, they exhibit a number of new effects due to the underlying non- 
uniform cllipticity of the family. The effective behaviour of such media in the vanishing period limit 
also differs notably from that of multi-dimensional models investigated thus far by other authors, 
due to the fact that neither component of the composite forms a connected set. We then discuss 
a modified problem, where the equation coefficient is set to a positive constant on an interval that 
is independent of the period. Formal asymptotic analysis and numerical tests with finite elements 
suggest the existence of localised eigenfunctions ( "defect modes" ) , whose eigenvalues situated in the 
gaps of the limit spectrum for the unperturbed problem. 
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1. Introduction. 

1.1. The general context for the problem in hand. The description of 
the effective behaviour of high-contrast composites ("high-contrast homogenisation") 
has been of particular interest in the analysis and applied communities over the last 
decade. The analytical part of the related literature starts with the work [T7], which 
developed in detail some earlier ideas of pQ concerning the use of "two-scale conver- 
gence" for the analysis of the limit behaviour of the boundary-value problem 

-div(A s (x/s)Vu) = f, /ei 2 (f>), ueH'(n), A e =e 2 Xo I + xiI, e > 0, 

where C K™ is a bounded domain, and \q, \i are the indicator functions of [0, l) n - 
periodic sets in R n such that Xo + Xi = 1- 

Several contributions to the high-contrast homogenisation followed: in the linear 
and non-linear, scalar and vector contexts, with various sets of assumptions about the 
underlying geometry of the composite. With applications mainly in solid mechanics 
and electromagnetism, high-contrast media have served as a theoretical ground for a 
number of effects observed in physics experiments, in particular those related to pho- 
tonic band-gap materials and cloaking metamaterials ([2]). The range of techniques 
developed in these contexts and their applications continue their rapid expansion, and 
the present paper is one contribution aimed at addressing some aspects that have thus 
far been left out of the scope of the related research. 

More specifically, we approach the question of the analysis of the spectral be- 
haviour of high-contrast composites in the case when the component represented by 
the function xi (the "matrix" of the composite) is disconnected in R™. Clearly, this 
is always the case in one dimension (n = 1), which is the situation we study in the 
present article. 
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1.2. Problem setup. We consider solutions u to the following family of elliptic 
problems on an interval (a, 6)cM: 

A e u-\u = f, feL 2 (a,b), £>0, AeC, (1.1) 

where the operators A e are given by the closed bilinear form 



(A e u,v)= p(x/e)(e 2 xo(x/e) + xi(x/e))u l (x)v'(x)dx } u,«6Jo. (1.2) 

J a 

Here p = p(y) > is a 1-periodic function in R such that € L°°(0, 1), the 

functions %o an d Xi are the indicator functions of 1-periodic open sets Fq and F± 
such that FqU Fi = R, and f) denotes a closed linear subspace of H 1 {a : b) that 
contains Cg°(a,b). We make no assumptions regarding boundedness of the interval 
(a, b), in particular it may coincide with the whole space R. 



In applied contexts the problem (1.1) corresponds to, e.g., the study of wave 



propagation in a layered 2D or 3D composite structure where / = 0, A > 0. In what 



follows we study the spectrum S £ of the problem ( 1.1 ), i.e. the set of values of A for 
which A 6 — XI does not have a bounded inverse in L 2 (a,b). Throughout the article 
we employ the notation a{A) for the spectrum of an operator A, and the notation 
Q for the "unit cell" [0, 1) whenever we describe the behaviour with respect to the 
"physical" variables x, y. We continue writing [0, 1) for the "Floquet-Bloch dual" cell 
when we refer to the domain of the quasimomentum 9. 

1.3. Our strategy for the analysis of ( |1.1[ ). It has been well understood 
in the existing literature on the subject (see [2]. |17j. [19j). that in the analysis of 
convergence of spectra of families of differential operators with periodic rapidly oscil- 
lating coefficient, one has to deal with two distinct issues: the lower semicontinuity 
of the spectra in the sense of Hausdorff convergence of sets, and the possibility of 
spectral pollution, the lack of which is often referred to as "spectral completeness". 
The former issue, which in the wider spectral analytic context has been looked at 
from a more general perspective (see e.g. [4]), is usually dealt with by proving first 
a variant of strong resolvent convergence. In the case of periodic operators involving 
multiple scales, one typically makes use of the so-called "two-scale convergence" (see 
e.g. [13], [T], [II])- In the present paper we follow this general approach in proving 
the related lower semicontinuity statements both for the whole-space problem and for 
the problem in a bounded interval. It should be pointed out that this first part of the 
analysis of spectral convergence is not completely independent from the subsequent 
study of spectral completeness: unless some assumptions are made concerning the 
geometry of the periodic composite in question (see e.g. [17] ) , one may not get the 
best possible "lower bound" for the limit spectrum. It has been noticed that, in order 
to capture the behaviour with respect to all Bloch components in the limit as e — > 0, 
it is preferable to use an advanced, "multi-cell" version, of the standard two-scale 
convergence; see e.g. [1], Chapter 5], where this more refined approach is adopted. 
It is a version of this last, more detailed, procedure that we adopt in the present 
article. 

In the proof of spectral completeness, a natural strategy seems to try and analyse 
the relative strength of different Bloch components in a given (convergent) sequence 
of eigenf unctions. This idea has been elaborated in [5| in the specific context of 
"high-frequency" homogenisation with the use of what the authors refer to as the 
"Bloch measures" . A combination of a compactness argument in the related space 
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of measures and a special "slow-variable modulation" construction then yields the 
simultaneous convergence of the given sequence to a limit eigenfunction and of the 
associated eigenvalues. In the present work we suggest an alternative approach (see 
Section [3]) to the convergence of eigenfunctions, which we believe is closer in spirit 
to the idea of "spectral compactness", i.e. compactness of eigenfunctions in a norm- 
preserving topology. Our approach is based on the idea that once one has control of 
the behaviour of eigenfunctions in the orthogonal complement to the space spanned 
by the limit eigenfunctions, one can immediately pass to the limit, as e — > 0, in the 
weak formulation of the original family of eigenvalue problems. This idea allows us 
to cover the analysis of spectral convergence for a wide range of operator families, 
including those considered by [2],[l7], [3 Chapter 4]. 

The key element in our analysis, which allows us to implement the above idea 



is Proposition 3.2 below (see Section [3.1[ ), or equivalently Proposition 3.3 These 
statements establish a uniform version of the Poincare-type inequality between the 
projection of a given function onto the "poorly behaving" subspace and the L 2 -norm 
of its derivative on the part of the domain where solutions of the eigenvalue problem 
can be shown to be a priori small as e -> 0. Different versions of the same idea have 
appeared in a number of other contexts, serving a similar purpose of "compensating" 
somehow the apparent loss of compactness in the problem, for example, in the form 
of Korn inequality in elasticity (see e.g. [BJ, and also [TS] for its multiscale versions), 
in the form of the so-called "energy method" in classical homogenisation (see [12] L 
and, more recently, in the form of a "generalised Weyl decomposition" for problems 
with degeneracies (see (HI])- For nonlinear variants of the same idea, the reader may 
be referred to the "geometric rigidity" (see [7j) and "^4-quasiconvexity" (see [5]). 
For an easier introduction to the problem, in what follows we start with the 



analysis of the problem (1.1 ) in the whole-space case, (a, b) = E, see Section [2j While 



a version of the the compactness argument developed in the bounded-interval setting 
(see Section |3| applies here as well (once complemented by a suitable Weyl-sequence 
argument), we present a different argument, based on some ideas of [SJ Chapter 5], 
where the spectral analysis is carried out in a more challenging setting of the Maxwell 
system. 

Throughout the article we assume for simplicity that the restriction of xo to the 
periodicity cell [0, 1) is the indicator function of an open interval (a, /3), which we also 
denote by Q . We use the notation Qi for the interior of the complement of Qo to 
the interval (0, 1). 

2. Limit analysis for the whole space. In this section we consider the case 
(a, b) = M. One well-known procedure for calculating S e is the Floquet-Bloch decom- 
position ([3]) following the rescaling y = x/e. Then, for 9 £ (0,1] the sequence of 
eigenvalues A = X(8) corresponding to 6>-quasiperiodic solutions to the Floquet-Bloch 
problem on the interval (0, 1) associated to the differential expression (p(e 2 Xo+Xi) u 'Y 
is obtained by solving the dispersion equation 



\ X(a - P + 1)) sin(VA(a - /?)) 



+ cos(eVA(a! + cos(yX(a~ P)) = cos(2tt6>). 
Passing to the limit in the above equation as e — > yields 

~(a-P+ l)VAsin(VX(a - /?)) + cos(VA(a - = cos(2tt(9). (2.1) 
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By varying 9 as indicated we obtain (for a = 1/4, f3 = 3/4) the set shown in Fig. 



2.1 




Fig. 2.1. The square root of the limit Block spectrum. The oscillating solid line is the graph of 
the function f(t) = cos(t/2) — t sin(t/2)/4, where t represents \f\ in the formula \2. l\ with a = 1/4, 
/3 = 3/4. The square root of the spectrum is the union of the intervals indicated by bold lines. 



Our first result is the following theorem. 

Theorem 2.1. Let (a, b) = R. Then the set lim^o S s is given by the union of 



solution sets for the equation (2.1) for all 9 G [0, 1). 

For 9 € [0, 1), we denote by Hg(Q) the space of functions u £ that are 

0-quasiperiodic, i.e. such that v(y) = exp(2iri9y)u(y) , y G Q, for somcQu e H#(Q)- 
We also denote 

V(9) := {v e Hg(Q) : p(y)v'(y) - for y e Qi} . (2.2) 
Consider an operator A(9) such that 



(A(9)u,<p)= P (y)u'(y)^(y)dy V<p£V(6), u G dom(A(0)) c V(6), 



defined on the maximal possible domain dom(A(9)). Henceforth (■, ■) denotes the usual 
inner product in L 2 (Q). By a standard argument (see e.g. |llj ) such an operator ex- 
ists, is unique and, under the adopted conditions on the coefficient p, is self-adjoint 
and has compact inverse (except for the case 9 = 0, when it has compact inverse as an 
operator on V{9) G C). Therefore the spectrum a(A(9)) is discrete and unbounded, 
i.e. it consists of eigenvalues < Ai(0) < \2(0) < of finite multiplicity with 
eigenfunctions v k {9) = v k (9,y). The eigenfunctions corresponding to different eigen- 
values are automatically orthogonal in L 2 (Q). We also carry out the orthogonalisation 
process on those eigenfunctions that correspond to the same eigenvalue, and normalise 
each eigenfunction so that ||^ fc (#)||L 2 (Q) = 1) f° r au ^ S [0, 1), k e N. 

Our aim is to show that the limit set lim £ _j.o S E coincides with the union of 
the spectra of the operators A(9), 9 £ [0,1), which, in turn, are described by the 



"dispersion relation" (2.1 1 



x As the notation Hq(Q) is usually reserved for the space of H X (Q) functions vanishing on the 
boundary of Q, we denote by H\(Q) the space Hl(Q) when 6 = 0. 
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In order to demonstrate first that the latter is included in the former, for each 
JV e N, we define an "intermediate" operator in L 2 (NQ), whose spectrum is 
contained in lim e _>.o S e and contains the spectrum of each of the operators A{&), 
6 = j/N, < j < N — 1. The details of this argument, which relies on a procedure 
that we refer to as "iVQ-periodic homogenisation" , are given in Appendix A. 

An essential component of the proof of the converse inclusion is the following 
lemma. 

Lemma 2.2. For any given 9 £ [0, 1) and ip £ V (9), let 9 £ £ [0, 1) be such that 
9 £ — > 9 as e — > 0. Then there exist p> £ £ V (9 e ) such that p> £ — » ip strongly in H (Q) 
as e — > 0. 

Proof. For 9 £ [0, 1) the space V{9) consists of functions that are #-quasiperiodic 
and constant in each connected component of Qx, that is for any ip £ V(9) one has 
f{v) = v) c + where c £ C, v £ Hy(a, b), and 

f 1, y£ [0,a), 

r)(0,y)~\ (exp(2m9)-l)(b-a)- 1 {y-a) + l, y£[a,b], 
[ exp(27ri0), y£ (6,1). 

For each value of e we now define ip e by the formula <p £ (y) = f](9e, y)c + v(y), y £ Q. 
Notice that by construction ip £ £ V (9 e ) and, since rj is uniformly continuous with 
respect to 9, one has ip 6 — > ip strongly in H 1 (Q). □ 

We next show that for a sequence A £ £ cr(A e ) such that X e —> A the inclusion 
A £ a {A (9)) holds for some 9 £ [0, 1). 

Theorem 2.3. Let \ e £ cr(A e ) such that X £ -> A. Then there exist 9 £ [0, 1) and 
u £ Hg(Q), u 7^ 0, such that 

( p(y)u'(yyp 7 (y)dy = X [ u(y)^y)dy, W^£V(9). (2.3) 

J a JO 

Proof. Since A e £ a(A £ ), by the Floquet-Bloch decomposition, there exists u £ £ 
Hg (e<9), u ^ 0, such that 



p(x/e)(e 2 xo(x/e) + xi{x / e)))u' e (x)<p' (x)dx = X £ / u £ (x)ip{x)dx (2.4) 

> JeQ 



for all tp £ Hi (eQ). Rescaling the formulation (2.4) with y = x/e yields the existence 



of u £ £ Hi (Q) 



' II u eIIl 2 (Q) 



= 1, such that 



e 2 / p(y)u' £ (y)(p'{y)dy + / p{y)u £ (y)ip{y)dy = X £ / u £ (y)ip(y)dy (2.5) 

Jqi Jq Jq 

for all ip £ Hl E (Q). 

The sequence 9 £ is bounded and therefore there exists some 9 £ [0, 1] such that, 
up to a subsequence which we do not relabel, 9 £ — > 9. Without loss of generality, if 
9 = 1 we set 9 — 0, so that 9 £ [0, 1). By substituting ip = u £ in ( ]2.5[ ), the sequence 
u £ satisfies the bounds 

II l?(q) - Cs > \\^ou' £ \\ L 2 {Q) <C (2.6) 

with a constant C > independent of e. 



() 
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Due to the weak compactness of bounded sets in H 1 {Q), the bounds (2.6), along 
with ||w e || i 2(Q) = 1, imply that, up to extracting a subsequence, u e converge weakly in 



H 1 (Q), and therefore strongly in L 2 (Q), to some u G H 1 (Q), 



\u \ 



l 2 (Q) 



1. Clearly, 



for w e (y) := exp(— 2iri9 £ y)u e {y) one has w £ G H^(Q), and the uniform convergence of 
exp(27ri6' e ?/) to exp(27ri#i/) as e — > implies that w E converge weakly in H 1 (Q) to wo 
given by the formula Wo(y) — exp(— 2-Ki6y)uo(y) , so that uq € Hg(Q). Furthermore, 
(2.6) implies that Xi u ' e ~^ strongly in L 2 (Q), hence uq G V (9) . 



In order to show that u satisfies the limit identity (2.3), for a fixed tpo G V (9) . 
let ip e G V (6 e ) be given by Lemma 2.2 Substituting ip e in (2.5), we obtain 



P{y)u' e {y)y' e {y)&y = A £ / u £ {y)ip £ (y)dy 



(2.7) 



By virtue of the facts that (p £ 



ipo strongly in H 1 (Q) and u £ — 1 u weakly in H 1 (Q), 
passing to the limit e — > in (2.7) immediately implies ( |2.3| . □ 

The above "limit spectrum" lim e _>o S e is strictly larger than the set obtained 
by the two-scale analysis of the operator A e of the paper [T7]. In particular, the 
spectrum of the homogenised operator obtained in [17] coincides with {Afc(0)}^ =1 , 
using our notation. Our analysis above shows that the set lim^o S e has, in fact, 
a band-gap structure, with infinitely many gaps opening in the interval [0, oo), as 
e — > 0. This fact suggests possible applications of the above composite structures 
to the design of optical or acoustic band-gap materials, which wc discuss in Section 
[4] The above effect also raises a mathematical question of the analysis of the limit 
behaviour of the operators A 6 in the case when (a, b) is a bounded interval, which we 
study in the next section. 

In what follows we assume that a, b G eFi. Our results are also easily carried over 



to the case when a, b G eFq, if one modifies (1.2) on those connected components of 
eFq that contain a or b, by changing the related coefficient from e 2 to unity. 

3. Spectral behaviour on a bounded interval. It is known that the classical, 



"moderate-contrast", analogue of the problem (1.1 )— ( 1.2 ) leads to limit spectra of 



different kinds for problems on bounded and unbounded intervals (a, b) : the limit set 
in the case of the problem in the whole space is purely absolutely continuous while 
in the case —00 < a < b < 00 it is purely discrete, i.e. it consists of eigenvalues 
with finite multiplicities, see e.g. [3]. A similar situation occurs in multidimensional 
high-contrast problems where the inclusion Fq D Q has a non-zero distance to the 
boundary of Q, see |17) . where, in addition, some eigenvalues of infinite multiplicity 
are present. 



As we shall see next, this is not the case for the problem ( 1.1 )-( 1.2 ), in particular 
rescaling y — x/e and replacing the form (1.2) with an integral over the whole of 
K, leads to higher-order errors in the limit as e — > 0, which can be ignored in the 
leading-order, "homogenised" , description of the operator A £ . 

In this section we employ, for convenience, the following notation: Q := (a,b), 

n e ■= n n (eFi), n e ■.= n n (eF ). 

3.1. The convergence result. The following theorem holds. 

Theorem 3.1. Consider an operator A e from the class described in Section \1J^ 
subject to the geometric modification mentioned at the end of Section [^[ The set 
lim^o S E is given by the union of solutions to the equation (2.1) for all 9 G [0, 1). In 



particular, it is independent of the choice of the space in Section 1.2 
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The inclusion of the union of the spectra cr(A(0)) in the set lim e _j.o S £ is proved 
in the same way as in the case of the whole-space problem, see the proof of Theorem 



2.1 In what follows we therefore discuss the converse inclusion (c/. Theorem 2.3). 



Notice first that for each A e € o-(A e ), there exists u e G fj, ||w e || i 2 



that 



(n) 



1, such 



p{^)u' e {x)ip'{x)dx + e 2 / p(f) u' E {x)y'(x)&x = X e / u e (x)<p(x)dx (3.1) 



for all if € f). Setting ip — u e in (3.1 ) yields the estimates ("a priori bounds" 



where C*b > is independent of e. 

For every bounded interval D we denote D\ 
space (c/. dpk, (53, d47}) 



ellz, 2 (ng) ^ ^b, (3.2) 
:= Z?nFi and introduce the function 



V(D) := {u e ff 1 ^) : fo) = /or y G 



The following statement is central to the proof of Theorem |3.1| 
Proposition 3.2. There exists a constant C± > such that 



P 



V(D)^u\\ Hl(D) < C±\\u'\\l*(Di) 

for any bounded interval D C K and any u € H 1 (D). Henceforth V(D) denotes 
the orthogonal complement of V(D) in the space H 1 {D) equipped with the usual 
inner product, and Py^±u denotes the orthogonal projection of the function u onto 
V(D)\ ^ 

Consider the norm 



u(y)dy 



1/2 



S(y)\ 2 dy 



(3.3) 



which, in view of Lemma 4.3 in Appendix B, is equivalent to the usual H -norm. 
In the same appendix the above Proposition |3.2| is shown to be equivalent to the 
following "uniform in 9" Poincare-type inequality. 

Proposition 3.3. There exists a constant C > 0, which depends on a and /3 
only, such that for any 9 g [0, 1) 



w 



< C 



\w'{y)\ 2 dy, Vwe^ffl). 



(3.4) 



We denote by V {&) the orthogonal complement oiV{9) defined by (2.2) in the space 
Hg(Q) equipped with the inner product in H 1 (Q) associated to the norm (3.3). 

Proof. We first note some properties of functions that belong to the space V ± (9), 
which follow immediately from the characterisation of the space V{9) given in the 



proof of Lemma 2.2 



Lemma 3.4. Let w € V ± (9), then 

(i) The equation w" {y) = holds for y € Qq. In particular, the function w is 

linear on the Qo-component of the unit cell: w(y) = (w(/3) — w{a))(j3 — a)^ 1 (y — 
w(a)) + w(a) for y G Q . 
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(ii) ForO^O one has (/? - a)- l (w(P) - w{a)) = (1 - cxp(-27ri6'))- 1 (a + (1 - 
p)ex.p(-2Trie))f Qi w(y)dy. 

(Hi) For 9 — one has Jq^ w(y)dy = 0. 

We now return to the proof of Proposition [33] We consider three different cases, 
depending on the location of the quasimomentum within the Floquet-Bloch cell [0, 1). 
Case I: 6 = 0. By (3.3) and Lemma [3.4| (i), (hi), we find that 



w 



\w'(y)\ 2 dy 



H/3) - w(a)\' 
P-a 



Since w(l) = w(0), we obtain the estimate 



w 



< / \w'(y)\ 2 dy + (\ w (p)- w (l)\ 2 + \w(0)-w(a)\ 2 ) 



< 1 



\w'(y)\ 2 dy 

, 2|Qi| 
P — a 



P-a 
2 

P-a 
\w'(y)\ 2 dy 



w'(y)dy 



w'{y)dy 



Case II: 9 e (0, S) U (1 — 6, 1), where < S < 1/2 is to be chosen appropriately. 
By (3.3) and Lemma 3.4 (i), (ii), we find that 



w(y)dy 



1 



f 



\w>(y)\ 2 dy+ lwm - w(a)] 

P — a 

\w(P) - w(a)\ 2 + [ \w'(y)fdy 



P-a \d g \< 

where d e := (P - a)(l - exp(-27ri6'))" 1 (Q; + (1 - p) exp(-2?ri6')). From the fact that 
w(l) = exp(27ri6')w(0) we infer 

\w(P) - w{a)\ 2 < 3 (HI) - w{P)\ 2 + \w{a) - w(0)\ 2 + \w(l) - w(0)\ : 



<3|Qi| / K(y)| 2 dy + 3|exp(2^) 
JQi 



i| 2 k(o)| 5 



and therefore 

IIMII 2 < (3|Qi| 



a 



\w'(y)\ dy 



+3 (/^ + wf) |expM - lfW0)12 

Notice that \d g \ 2 = (/3-a) 2 (2-2cos(27r6»))- 1 (a 2 + (l-/3) 2 + 2a(l-/3) cos(27r6»)), hence 
\ds\ vanishes at 6 — 1/2 for the special case a = 1 — p. In view of this observation 
and in order to have a bound on the constant de we require that 6 < 1/4. Further, 
by continuity of the embedding of H 1 (Q) in C(Q), there exists a constant c, which is 
independent of 9, such that 

ko(0)| < c|||tu|||, 
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and thus 



Mil 2 < (^1 + ^ + ^)/ gi k(,)l 2 d y 



1 J « . . . - ,2 -2||| 1 1 1 2 



H^+wJ 1 ^ 2 "^- 11 s|IHr (3 - 5) 

We now choose <5 < 1/4 so that ((£ - a)" 1 + \d e \~ 2 ) |exp(2 7 ri6') - 1| 2 c 2 < 1/2, and 



hence |(ie| _2 is bounded above by a constant independent of 0. The inequality (3.5 1 
now immediately implies the required estimate. 

Case III: 9 £ [6,1 — 8], For given x € (/?, 1], y € [0, a) we write 

c-x pa 

w(x) = v/(t)dt + w(P), w (y) = - w' {t)dt + w(a) 



which implies, in view of Proposition 3.4 (ii), 



w(x) - w(y) = w(P) - w(a) + ( J +J ) w'(t)dt 



X />Q 



= d e / w(y)dy+ / + / u/(t)df. 

JQi V.//3 Jy J 

In particular, substituting x = 1, y = and using the fact that w{l) = exp(27ri^)w(0), 
we obtain 

W W = 1 T-V^m / w (y) d y+l ! / w'(y)dy, 

1 - cxp(-27ri6 l ) y Ql 1 - exp(-27ri6>) J Ql 

whence 

w(x) = - [ w'(t)dt + w{l) = - [ w'(t)dt 



d f 1 

1 — exp(— 27T1&J Jq 1 1 — cxp(— 27rit/ 



Qi 



Integrating the last identity over (/3, 1] yields 
w(y)dy = - [ ( [ w'(t)dt] dx 



p J/3 \Jx 



1 — exp(— 27ri&) 1 — exp(— 27ri0 



Qi 



Similarly, we write 
rv 

w ( y )= w ' (t)dt + exp(-27ri0)ti;(l) 
./o 



/■» , exp(-27rig)d 9 /" exp(-27rig) f 

= / w {t)dt + / w(y dy + / to (y)dy, 

Jo 1 - cxp(-27ri0) 7 Qi 1 - exp(-27ri6») J Ql 
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which upon integration over (0, a) yields 
w(y)dy = ^ U\'(t)d?jdy 



aexp(—2ni9)dg f , , , aexp(— 2ni9) f ,. 

/ w{y)dy+- / uj(j/)dy. 

1 - exp(-27ri6>) y Q 1 - exp(-27ri0) y Q 



(3.7) 



Combining equations (3.6) and (3.7) we obtain 
1 - — — ' " ~~7\ — 1 / w(y)dy 



1 - exp(-27ri6l) 



a / ry 



WO 



u/(i)dt dy 



pKJx'j i - cxp(-27ri0) y Ql w 

Squaring both sides and using the Cauchy-Schwarz inequality yields 



1 - 



(1 - p + a exp(-27ri0)) dg 



1 — exp(— 27ri6 



w(y)dy 



< 2 4 



|l-/3 + aexp(-27ri6»)| 2 

|l-cxp(-2m0)|2 M Q 



K(2/)| 2 dy. 



A direct calculation shows that the coefficient in the left-hand side of the last inequal- 
ity is separated from zero in the range of 9 considered. 
Finally, we argue that 

\w(fl) - w(a)\ < 4(H/3) - iu(l)| + \w(0) - w(a)\ + \w(0)\ + \w(l)\) 



and 



w'{y)dy 



w'(y)dy 



+ \ W (0)\ + \ W (1)\ 



(l-/9)u;(l)= / / w'(t)dtdx+ [ w(x)dx, 

J p Jx J p 



aw(0) 



w'(t)dtdx + / w(x)dx. 
Jo Jo 



The required inequality follows, since by ( |3.3| and Lemma 3.4 (i), 

2 



Qi 



w(y)dy 



\w(P) - w{a)Y 



\w'(y)\ 2 dy. 



This completes the proof of Proposition |3.3| □ 
We now resume the proof of Theorem |3.1| 
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Let N e be the smallest integer such that O C e[a/e] + U^lQ 1 e(n,n + 1) =: fi e . 
For each e, we denote D e := e _1 f2 e and use the rescaling operator T e defined by 
(T e u)(y) = u(ey), y G D £ , and an extension operator E e from f) to H 1 (Cl e ) such that 

\\E(u £ )\\ L2(M < (1 + C7 E e)KIU»(n) , \\E(u e y\\ LH fr) < (1 + C E e)|K|U a( n) (3.8) 

for all e > 0, where the constant Ce > is the same for all e. (It is clear, for example, 
that the continuous extension of u e by constants satisfies this requirement.) We then 
define functions U e by the formula U e := T" 1 P v me\T s E(u s ) for each value of e. 

Extending each of the functions U s to the whole of M with period eN £ , we consider 
the "discrete Floquet-Bloch transform" of the functions U £ as follows 

u tW = aT E Ue{x + ek) e M 2 ^jk/N e ), x g n e , j = o, i, 7v £ - i, 

6 fc=0 

so that T e £/| G V(j/N e ) and U E = X^o^ ^1- Using the fact that the eigenfunctions 
v k (j/N e ), k G N, form a complete system in the L 2 (Q)-closure of the set V(j/N e ), we 
represent C/| in terms of them so that 

N c -1 oo 

Ue(x)= E IWlfcKGfc.!). (3.9) 

j=0 k=l 



N e -l oo 

c/^) = E - 1 e E^(it)(- fc )'( 

j=0 fe=l 



J «5 



(3.10) 



where U k (j/N e ) G C. Applying the Parseval identity to (3.9), followed by Proposition 



3.2 and the first of the a priori bounds (3.2) yields 

JV.-l oo 



j=0 k = l 



The last estimate, in combination with the first inequality in (3.8) and the fact that 
IK || z,2 (n) = 1, implies 



N E -l oo 

l-C ± C B e<eN e E E|^( 

3=0 k=l 



< 1 + C E £ 



(3.H) 



Denoting by <5(- — 9) the Dirac mass at 8, we infer from (3.11 1 the existence of C > 
such that 



fc=i 



N E -1 



E / d ^ - 1 < ^, d/x*(0) := eN E E K 9 ~±)M- ( 3 - 12 ) 



Clearly, for each fc the sequence {/^je is bounded in the space of Radon measures on 
[0, 1]. Therefore, up to a subsequence, fi k weakly converge as e — > to some measure 



[ ud^-> [ udfi k , Vu6C[0,l] 
Jo Jo 



(3.13) 
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The above result follows from recalling that the space of finite Radon measures on 
[0, 1] coincides with the dual space C[0, 1]* and hence bounded sets of Radon measures 
and relatively compact with respect to weak star convergence in this space. 



the second of the a priori bounds (3.2), we write 



Furthermore, taking into account (3.10), the second of the estimates (3.8) and 



N e -\ oo 



c B (i+c B£ )> *m\>ito=*N. E E|^(*)| IK«*)'Gfc.O 



3=0 k=l 



L2 (Q) 



(3.14) 



N s -1 

> \\ P -i\\ L ~ iQ) sN e E|^(i) 

j=0 fe=l 



(3.15) 



oo -i 

= lb- 1 IU~ ( o)E / W)*n k M- 
k=l J ° 



(3.16) 



In order to obtain (3.15), we use the fact that v k (9), k e N, are eigenfunctions 
of the operator A(9) with eigenvalues Xk(0), and (3.16) is a version of the same 
expression written in terms of the measures defined in (3.12). For any integer K > 2 



the inequality ( 3.14 )-( 3.16) immediately implies 



k=K 



£ / ^e(6)<\\p-%L Q) C B (l + C E e) mm \ K (0) 



-I i\— >oo 



0. 



In view of (3.12) and (3.13) we thus argue that Y^k=i M fc ([0; 1]) = 1) hence fi k is a 
non-zero measure for some k. 

We next show that A = Afc (#) for some ko, 6. To this end we set ip = U e in the 
weak formulation (3.1): 



p (^) u' e (x)U^(x)dx + e 2 p (§) u' e {x)U' e (x)dx = X e u e {x)U e (x)dx. (3.17) 

Of J!H is] 



Notice that, by the definition of the space V{D e ), the first term in the equation (3.17) 
vanishes. We also have, in view of Proposition |3.2| and the second of the estimates 



(3.2), as e -> 



-/ p(^)u' e {x)U' £ {x)dx = e 2 / p(f)U' £ {x)U' s {x)dx + o{l) 



oo oo N E — 1 

EE E 

fc=l /=l j=0 "' £ yo 



p(f)^(]t)( w *)'(it.f) r3 2(it)( w ')'(ifc'f) da! + ( 1 ) 



AT e — 1 



^EE A ^(i)F(i)| +°( 1 ) = E/ A*(e)d^ + (i) 

fc=l 1=0 fe=l ^° 



(3.18) 
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In the last three equalities of (3.18) we use (3.9), ( 3.12[ ) and the fact that v k {9) are 
orthogonal eigenfunctions of A(8) with eigenvalues Xk(9). Similarly we have 



/ u e (x)U £ (x)dx = Vd/4 + o(l), 
70 fe=i 



as e 0. Combining (3.17), (3.181, (3.19) yields 



oo „i oo -X 

2 / A fe (0)d/4(0)=A £ £ / d/x* 



o(l). 



(3.19) 



(3.20) 



Finally, passing to the limit e — >• in (3.20) and using the fact A& are continuous 
functions of 9 (see Appendix C) yields 



oo „i oo «x 

2 / A fc (0)dM fc (0) = A£ / 
k=i Jo k=i Jo 



d// 



We already know that /i fe ° is a non-zero measure for some fco, hence A = Xk a (0) for 
some 9. 

4. A modified problem with a compact perturbation, and the associ- 
ated defect modes. 

4.1. Ana lytical setup. In this section we discuss a modified version of the setup 
of Section 



1.2 



as follows. Consider the operator A £ defined by the bilinear form 

b 



(A e u,v) = j (pdXd(x) + p{x/e)(e 2 xo(x/e) + xi(x/e)) (l - Xd(x))ju' (x)v' (x)dx, 

(4.1) 

u,v G $j. Here the space Sj is as before, and Xd is the indicator function of a "defect" 
interval Id whose closure is assumed to be contained in (a, b), and pd is the corre- 
sponding "defect" coefficient (or "defect strength"). Analogously to the above, we 
assume that a, b, and the end-points of Ja belong to the set eF\, otherwise we modify 
(4.1) on those connected components of eFq that contain a, b and the end-points of 
Id by changing the related coefficient from e 2 to the unity. We denote by S e the 
spectrum of the operator A e . 

A formal two-scale asymptotic procedure carried out on the equation (cf. (1.1)) 
A e u = \u suggests that: 

1) The set lim e _j.o S £ is independent of the choice of the space fj and is given by 
the union of solutions to the equation (2.1 ) for all 9 £ [0, 1) and a sequence of "defect 
eigenvalues" {pd7r 2 j 2 /|^d| 2 }°° 



j=i> 

2) The "defect eigenfunctions" Uj corresponding to the above eigenvalues, j € N, 
decay exponentially away from the boundary of the defect: 

\uj(x)\ < Cexp(-e _1 dist{a;, J d }), x€(a,b)\I d , C > 0. 

We next present numerical evidence that supports these claims. 

4.2. Numerical results for the modified problem. We consider a defect of 
length |Jd| = 1/2 and strength p d = 2 in the middle of the interval (a,b) = (0, 1), 
i.e. Id = (1/4,3/4). For each value of e such that N := is a positive integer, we 
describe the intervals (0, 1/4) and (3/4, 1) on either sides of the defect according to 
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the equation (4.1 1: each of them consists of N cells of the same length l/(4iV), and 
in one half of each cell the coefficient in the form (the "modulus") takes the value 
1/(AN) 2 while in the other half it is equal to unity. Importantly, we assume that 
the endpoints both of the interval (a, 6) and of the defect are belong to the material 
component where the modulus is equal to unity. 

The results of solving the above problems with finite elements are given in Tables 
4.1 and |4.2| The values for the trapped mode are in good agreement with the values 



obtained by the asymptotic method: Ai 2) = 78.9568, \l°> = 315.8273, Ai° ; = 710.6115, 
Ai 6 ' = 1263.3094, where the superscript is the number of the stop band containing 
the related eigenvalue. 

In addition, the profiles obtained for such trapped modes (see Figure 4.1 for the 
case of periodic boundary conditions) suggest that the number of half-oscillations in 
a trapped mode is equal to the number of the mode in the sequence, which resembles 
the behaviour of the usual Neumann eigenfunctions on the defect. We also note that 
the decay of the trapped modes appears to be exponential, as can be seen in Figure 
4.2 the larger the contrast (and hence the number of subdivisions of the string) the 
more localised the mode, irrespective of the boundary conditions at the endpoints of 
the string. 



(3) _ 



k(«0 - 



Dirichlet boundary conditions 


, (fc,128) 
min 


. (fc,128) 
^max 


A* 


, (fc,256) 
A* 


. (fc,512) 
A* 


, (fc,1024) 
A* 


11.7939 


39.4603 










65.7875 


157.8859 


75.7674 


77.2502 


78.0741 


78.7304 


187.6799 


355.2599 


293.9534 


304.1141 


309.7163 


314.2461 


386.1413 


622.2747 










662.9213 


986.7685 


682.6577 


694.4984 


702.0486 


708.4576 


1018.4394 


1421.0468 


1225.1298 


1232.2190 


1243.1182 


1258.2799 



Table 4.1 

Stop bands and trapped modes for the modified problem with a defect, subject to the Dirichlet 
boundary conditions: A' fc ' 128 ' an d A^x 28 '' are the lower and upper bounds of the k th stop band for 



N = 128, and X\ 



A 



(fe,512) ,(fc,1024) 



are the trapped-mode eigenvalues in the k th stop 



band, evaluated for N = 128, N = 256, N = 512 and N = 1024 respectively. 



Neumann boundary conditions 


, (fc,128) 
min 


. (fe,128) 
A max 


, (fc,128) 
A* 


, (fe,256) 
A* 


, (fe.512) 
A* 


> (fc,1024) 
A* 


11.7515 


39.4980 










65.8359 


157.8901 


75.7676 


77.2509 


78.0741 


78.7314 


187.7334 


355.2765 


293.9539 


304.1145 


309.7164 


314.3057 


386.2091 


622.2747 










662.9779 


986.7698 


682.6578 


694.4985 


702.0486 


708.6496 


1018.5163 


1421.0556 


1225.1298 


1232.2190 


1243.1193 


1258.2626 



Table 4.2 



Stop bands and trapped modes for the modified problem with a defect, subject to the Neumann 
boundary conditions: A^j'^ 28 ^ and A^ix 28 ^ ore the lower and upper bounds of the k th stop band for 
N = 128, and \+ k,12S \ A^' 256 ', A^ fe ' 512 ', A^' 1024 ' are the trapped-mode eigenvalues in the k th stop 
band, evaluated for N = 128, N = 256, N = 512 and N = 1024 respectively. 
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(a) 



(fa) 




Fig. 4.1. Trapped eigenmodes for the modified problem with a defect and periodic boundary 



conditions (N=1024), corresponding to eigenfrequencies: (a) A 
314.24, (c) a1 5 ' 1024) = 708.46, (d) a1 6 ' 1024) = 1258.28. 



(2,1024) 



78.07, (b) A 



(3,1024) 



(a) 



(b) 




Fig. 4.2. Decay of the third trapped eigenmode (located in the fifth stop band) for the modified 
problem with a defect and periodic boundary conditions, as a function of contrast: (a) A^ 5 ' 32 ' = 



668.86, (b) A 



(5,64) 



682.89. 



(c) Ai 5 < 128) = 694.50, (d) A™ 



702.68. 



4.3. Photonic band gaps and trapped modes in high-contrast multi- 
layered dielectric structures. The string problem emerges, among other contexts, 
in the study of wave propagation in one-dimensional photonic crystals i.e. multi- 
layered dielectric structures invariant along two directions. In what follows we set 
these directions to be X\ and X3 in the usual Euclidean representation x = (x±, X2, £3). 
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We consider those solutions (£,%) to the classical system of Maxwell equations 
(0) that have the form 

£(x 1 ,x 2 , X3, t) — E(x 2 ) exp(z(ra 3 — ujt)), H(xi,x 2 ,x 3 , t) = H{x 2 ) exp(z(ra 3 — ujt)), 

where t is time, uj is the angular frequency, and n > is a "propagation constant". 
We write the Maxwell equations for the field variable {E, H) : 

E' 3 — inE 2 =iuifj,Hi, [ — H'3 + ikH 2 = iwsEi, 
ikEi = \bj[iH 2l < —inHi = iujeE 2 , 



-E[ = icj/iiJ 3 , I H[ = iueE; 



3- 



Here (i is the magnetic permeability, e is the electric permittivity at each point of the 
dielectric. 

We rearrange the above six equations in two groups of equations for (£1, H 2 , H 3 ) 
(transverse magnetic polarisation), and {Hi, E 2 ,E 3 ) (transverse electric polarisation). 
We choose Ei and Hi as the unknown functions within the respective groups and no- 
tice that the remaining unknowns are expressed in terms of these two scalar functions 
only. The equations satisfied by E\, Hi are 

{E[)' + {u 2 ne-K 2 )Ei =0, (4.2) 

(^ffD' + ^/i-r^ji/^O. (4.3) 
Note that (4.3 1 coincides with ( |1.1[ ) when n = 0, by setting 



oj 2 = A, (i = l, e 1 {x 2 )=p{x 2 /v)(v 2 Xo{x2/v)+Xi{x 2 /v)), x 2 E {a,b), i] > 0, 

where we use rj rather than e to denote the structure period, in order to avoid confusion 
with the standard notation for electric permittivity. Our analysis in Sections [2] and [3] 
carries over to the case k > 0, where we get a K-dependent version of the dispersion 



relation (2.1 1, as follows: 

i(a - $ + 1) - sin(y\{a - /?)) + cos(V\{a - /?)) = cos(2vr6i). 

Assuming infinitely conducting walls on either side of the dielectric (see e.g. [20] for 
further details), we supply (4.2) and (4.3) with homogeneous Dirichlet and Neumann 
boundary conditions respectively. 

In the numerical solution of above problem we employ finite elements with per- 
fectly matched layers, i.e. anisotropic absorptive reflectionless layers (see e.g. [20 ), 
on the top and bottom of the computational domain. Our results are shown in Figure 
|4.3| for k = .1, N — 16, and for the transverse electric mode with frequency A* = 78.34 
inside the third stop band. The latter corresponds to the first trapped mode shown in 



Figure 4.1 'a), in view of the fact that for k = .1 there is an additional zero- frequency 
stop band. The magnetic component of this mode (Figures |4.3[ a) and (d)) clearly 
shares the same features as the string mode in Figure |4.1[a) . 



Appendix A. In this appendix we argue that for any fixed TV G N and A < 0, 
say A = — 1, the solutions to the problems ( |1.1[ ) converge in an appropriate two- 
scale sense (see e.g. [1]) to the solution of some limit problem parametrised by N. 
Throughout our argument, we assume that (a, b) = R, although the results are valid 
for an arbitrary interval (a, b) C R, irrespective of the boundary conditions. 

We first formulate the related statement for N = 1. 
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Fig. 4.3. Obliquely propagating transverse electric wave in a high-contrast dielectric multilay- 
ered planar waveguide with infinitely conducting walls: (a) 2D plot of Hi; (b) 2D plot of E2; (c) 
2D plot of E3; (d) Profile of Hi along the horizontal centreline. Here N = 16, k = .1 and the 
normalized frequency A* = 78.34. 



4.4. Periodic homogenisation. Lemma 4.1. Set A = — 1 and let u e g H X {R) 
be the solution to the problem (1.1). Then there exists u(x,y) £ L 2 (fl; Vi) such that 



u(x,y), 



eu £ u v (x,y), 



xi(!M!K(*)^o. 



Here the space V\ is defined by 

Vi := {v E Hl{Q) : p(y)v'(y) = for y g Qi} . 
The function u G L 2 (R; V\) is the unique solution to the problem 



p{y)u yV (x,y)ip'(y)ip(x) dydx 



u(x,y)ip(y)ip(x) dydx 



(4.4) 



f(xMy)^(x) dydx, Vcp g C °°(M), ^ € V x . 



(4.5) 



Q JQ 



Note that the x dependence in the formulation (4.5| is trivial, indeed u(x,y) 



f(x)v(y) where v(y) G V\ is the unique solution to the problem 



Qn 



PivW (y)(fi' {y)dy + / v(y)ip(y)dy = / (fi(y)dy, Vip g V x . 



(4.6) 



This observation implies, in particular, that the spectra corresponding to (4.5) and 
(4.6) coincide. Therefore, denoting by A\ the operator defined by the form 

bi(u,v)= p(y)u'(y)v'(y)dy, u,v&Vi, 



we obtain lim e _>.o S e D a (Ai) . 
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4.5. 'WQ-periodic" homogenisation. In the argument above we found the 
two-scale limit operator A\ by choosing the periodic reference cell to be Q and passing 



analogue of Lemma 4.1 as follows 



to the two-scale limit in ( 1.1 1 as e — >• 0. Replacing Q with NQ, N G N, we obtain an 



Lemma 4.2. Set A = —1 and let u £ be the solution to (1.1). Then u £ — v un , 
up to some sequence we do not relabel, where un = f(x)vN(y) and vn is the unique 
solution to 

Xo(y)p(y)(v N )'(y)<p'(y)dy + v N (y)(p(y)dy = / f(y)dy, V<p G Viv- 
nq Jnq Jnq 

Here we denote 

V N := {v g H^NQ) : p(y)v'(y) = for y g Q,} . (4.7) 

Furthermore, lim £ _>o S s D a (An) , where An is the operator defined using the bilinear 
form 



b N (u,v):= Xo(y)p(y)u'(y)v'(y)dy, u,v € H#(NQ). (4.8) 

JNQ 



Applying Lemma 4.2 for all N g N yields 

]ima(A s ) D M a(A N ). 

NeN 

4.6. Relation to the Bloch spectrum. Notice that if 9 = j/N for some inte- 
gers N,j such that < j < N — 1, then all eigenfunctions v k (9), k = 1, 2, are iV- 
periodic, in particular v k (9) g Vjv. In fact, v k (9), k — 1, 2, ... are eigenfunctions of Ajv. 
Indeed, for any fixed f G Vn define the function $(y) := ^fc^o 1 viv + k) exp(— 2iri9k) 
and notice that since g one has (A(6»)w fc (6»), $) = A fc (0)(t>*(0), There- 

fore, writing for brevity v k (9) — v, we obtain 

JV-i . 

Xo(y)p(y)v'(y)cp(y)dy=y2 / xo(y + k)p(y + k)v'(y + k)ip(y + k)dy 

NQ k=Q JQ 



Xo(y)p(y)v(y)${y)dy = A fe (0) / v(y)<l>(y)dy 



Q 



w -! r r 

= >"Afc(0) / exp(27ri^)^(y + k)dy = A fc (0) / v(y)<p(y)dy. 

The above observations show that 

limS e D |J a {An) D \J a(A(j/N)). 

NeN New 

0<j<N-l 

Using the facts that the set of rational numbers j/N is dense in [0, 1) and that the 
eigenvalues A = X(9) are continuous with respect to 9 (see Appendix C below) yields 

hm^D |J a(A(9)). 

66 [0,1) 
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Appendix B. 

4.7. One classical inequality. Here, for reader's convenience, we give the proof 
of a version of the classical Poincare inequality (see e.g. [B]), which we use in the 
present paper. 

Lemma 4.3. There exists a positive constant Cp, which depends on a and ft only, 
such that 



Hy)\ 2 dy<C P 



u(y)dy 



+ J^\u'(y)\ 2 dy\ , Mu&H l {Q). (4.9) 



Proof. For fixed x € Q, y € Q%, we have \x — y\ < 1 and 

u(x) — u{y) = I u'(t)dt, 

which implies 

\u{x)\ 2 + \u{y)\ 2 ~2^(u{x)^y)) - 
Integrating hrst with respect to x, then with respect to y, yields 
\Qi\ [ \u(x)\ 2 dx+ [ \u(y)\ 2 dy 



u'(t)dt 



< \x-y\ I \u'{t)\ z dt< 
y 



u'(y)\ 2 dy. 



< 23? 



u(x)dx / u(y)dy 



+ \Q 1 \ / \u'{y)\ 2 dy. 



Here we have used the fact that \Q\ = 1. Now using the inequalities (4.7), ab < 
ea 2 + b 2 /e for any real a,b, e > 0, and 



we obtain 



IQil / l«(*)l 

Jo 



2dx+ w\ 



u(y)dy 



u(y)dy 



< \Q ± \ / \u(y)\ 2 dy, 



:( / l«(y)| a d»J+g 



< |Q X | / \u(x)\ 2 dx+ I \u(y)\ 2 dy 
Q 

2 



u(y)dy 



H-lQil / \u'(y)\ 2 dy. 



Setting e = \Q\\/2 gives 



\Qi 



^ I K^)| 2 dx<-i 



IQi 



u(y)dy 



+ / \u'(y)\ 2 dy. 



This is (49} for C P = where |Qi| = 1 - p + a. D 
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4.8. The equivalence of Proposition |3.2| and Proposition |3.3[ Suppose 
that Proposition |3.3| holds and consi der a bounded interval D C K and a function 
u £ V(D) . For the proof of Theorem 3.1 it is sufficient to consider the case when D 



is an interval of integer length. However, for completeness we carry out the argument 
for an arbitrary D. 

Clearly, there is a positive integer N and an interval I := (I, I + N) D D, I £ M 
such that {1, 1 + 1} cfj. We extend the function u to a function that is periodic on 
the interval / and is such that u £ V(I) ± , keeping the same notation u for such an 
extension. 

Notice that for each j = 0, 1, N — 1 the function 

N 



'(x) = N^^uix-l + tycxpi^mjk/N) 



k=l 

belongs to the space V(9j) ± , 9j := j/N. Indeed, for any v £ V(6j) C Hg.(Q) one has 

v(x) — v#(x) exp(i6jx), (4-10) 
where u# £ H#(Q). We extend v# by periodicity to the whole of K and also extend v 



to the whole of R so that the formula (4.10 ) holds for any Then, for the extended 



function v, one clearly has v(- + I) £ V(I), and v(x — k) = v(x) exp(— 2-Kijk/N) for 
any i£i Therefore, for any j = 0, 1, N — 1, one has 

, N-l .i 

/ u 3 (x)v(x)dx = N^ 1 u(x~l + k) exp(-2vrijfc/iV) 



N-l r k+l 



v(x)dx 



*' / rK+x 

N' 1 ^ / u(x)exp(-2-Kijk/N)v(x - k)da 

I n Jk 



= N^ 1 J u(x - l)v(x)dx = N^ 1 J u(x)v(x + l)dx = 0. 
Now, using the Parseval identity and Lemma |4.3[ we obtain 

N N 

i=i j=x 



N 

< (C P + l)C 2 iV^ ||K)'||| 2(Qi) = (Cp + l)C 2 h'||i 2(7i) , 



where I\ := I D F±. By ensuring, when performing the above extension from D to J, 
that the inequality H^'Hfg/j^ < C|I u 'I!l 2 (Di) ^o\ds with a positive constant C that 
does not depend on the interval D 1 we complete the proof of the implication. 



Conversely, assume the validity of Proposition 3.2 and suppose that 9 £ [0, 1) is 



rational, N9 £ N. For any w £ V(9) , we notice that 



w(x) = u(x) exp(i9x) 



(4.11) 
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for some u G Hg(Q), and extend first the function u by periodicity to the interval 



D = NQ and then the function w according to the formula (4.11 1 to the same interval 
D. Then using an argument similar to the above it is shown that w G V(D) , and 
hence 

IM| H i(£>) < C|k'IU 2 (-Di) 



for some positive constant C. This automatically implies (3.4) since the norm of w is 
the same for any interval of length one. By continuity in 9 the statement is extended 
to establish the existence of C > that serves all 9 G [0, 1). 

5. Appendix C. The continuity of the family V(9) implies that for a given k 
dimensional subspace of V{9) we can find, for 9' close to 9, a k dimensional subspace 
of V{9') close to V(9), in the sense of [H]. More precisely, 

PROPOSITION 5.1. Let 9i e [0, 1) and Fi C V(9), dimFi = k. For all e > there 
exists S > such that for all 92, 1 6*2 — #i| < 5, there exists F2 C V{92), dimi 7 ^ = k 
such that max{dist(i 7 'i, F%), dist(i 7 2, Fi)} < £, where the distance between two linear 
subspaces N, M, of H 1 (Q), is defined by the formula 

dist(iV,M):= sup inf||u-i>||. 
ueN »e M 
NI H i=i 

Henceforth within this appendix, we write H 1 instead of H l (Q) for brevity. 

Proof. Let F\ C V(9i), dimi^ — k, and let n = 1, . . . , k, be a basis of Fx. 



Lemma 2.2 implies that for any 6*2 that is sufficiently close to 9\, there exist v n G V{92) 



that are close to in the iF-norm. We construct an orthonormal sequence as follows: 

n-1 

fl = V\, fi = 021 /1 + «2, = ^2 a-rimfm + v n, fl = 3, . . . , k, 

m=l 

where a nrn = — ("m /rrJ-ff 1 ■ By construction, fl is close to fl. Notice further 

that a,2i — \\vi\\ H \ {v2,vi)h 1 can be made as small as necessary if 92 is sufficiently 
close to 6*i . This implies that /f is close to f\ . By induction we show that is close 
to for all n = 1, 2, fc. 

Defining F 2 = span {/^ : n = 1, . . . , A;}, we argue that F 2 satisfies the properties 
of the proposition. Indeed, by construction, F% C ^(#2), dimF 2 = k. Further, for fixed 
u G Fi one has u = J2n=i ^n/ni f° r some 6„ G C. The function d = J2n=i bnfn belongs 
to F2 and is close to u if 6*2 is sufficiently close to 6\. It follows that dist(Fi, F2) < e 
when \&2 — 9i \ < S for sufficiently small 6. Reversing the roles of u and v completes 
the proof. □ 

Theorem 5.2. The eigenvalues Afc(0), k G N are continuous functions of 9 € 
[0,1). 

Proof. By a variational argument it is known that (see e.g. |15j ) 

A fe (0 = inf sup \. V ||2 ; ; , 5.1 
fcy(e) u6F \\u\l 2(0) 

dimF=fc W; 

Let Fi C V{9i), dimFi = fc. For a fixed e > 0, let S, 9 2 and F 2 C y(6> 2 ) be given by 
Proposition 5.1 For a fixed /2 G F2, there exists /1 G Fi such that ||/i — /2II//1 < £• 
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This implies 



(A(9 2 )f 2 J 2 ) (AjfiMUi 



\\M\lhq) II/iIIl2(q) 
It follows, by the arbitrary choice of f 2 , that 



{A{e x )u,u) i 

e < sup j 1" e - 

ue-Fi |M| i2(Q) 



sup 

«SF 2 



(^2) 



//. .7 ] (A(6i)u, u) 

< SUP ^--72 " 



E. 



i 2 (Q) 



L 2 (Q) 



Therefore, by the fact that F\ is arbitrary and in view of the equation (5.1 1, we obtain 



{A{9 2 )u,u) . 

mf sup ^ < mi sup 

F>«BF a \\U\\ 2 L2(Q) FlCVpiW, 

V ^ ! dim_Fi— k 



+ e 



li 2 (Q) 



= \ k (0)+s. (5.2) 
The set {F 2 } is a subset of the set of all k dimensional subspaces of V(9 2 ), therefore 



Xk(0 2 ) 



(A{6)u,u) (A(0 2 )u,u) 
vox sup ^ < mi sup 



ry ^ 1111 ULIJJ r\ 

FcV(6 2 ) ue F ||uH r2 ^ F2 „ e F 2 

dimF=fc L 



(5.3) 



L^(Q) 



Equations (5.2) and (5.3) imply 



\k(e 2 )-\k{0i)<e. 

Reversing the roles of d\ and 8 2 yields the desired result. □ 

The above statement of continuity of = \k(0), k £ N, is not a simple con- 
sequence of the continuity of eigenvalues for the usual Floquet-Bloch decomposition. 
An important distinct feature of the present statement is the dependence on of the 
operator domain V(6) in the variational principle (5.1). 
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